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, Abstract. Using two simple and robust inputs to constrain tlie mixing matrix of the isosinglet scalar mesons /()(1710), 

^ ■ /o(1500), /o(1370), we have shown that in the SU(3) symmetry limit, /o(1500) becomes a pure SU(3) octet and is degenerate 

' with ao(1450), while /o(1370) is mainly an SU(3) singlet with a slight mixing with the scalar glueball which is the primary 

^ I component of /o(1710). These features remain essentially unchanged even when SU(3) breaking is taken into account. We 

1^ , have deduced the mass of the pseudoscalar glueball G from an T]-T]'-G mixing formalism based on the anomalous Ward 

■ identity for transition matrix elements. With the inputs from the recent KLOE experiment, we found a solution for the 

' pseudoscalar glueball mass around (1.4±0.1) GeV. This affirms that T)(1405), having a large production rate in the radiative 

0^ I i/)// decay and not seen in yy reactions, is indeed a leading candidate for the pseudoscalar glueball. It is much lower than the 

results from quenched lattice QCD (> 2.0 GeV). 
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INTRODUCTION 



The existence of glueballs is a unique prediction of QCD as a confining theory. The latest lattice calculation of the 
glueball spectroscopy shows that the lightest glueballs are scalar, tensor and pseudoscalar glueballs with masses of 
order 1710, 2390 and 2560 MeV, respectively 1 1]. Since the lattice calculation was done in the quenched approxima- 
tion, the predicted masses are for pure glueballs in the Yang-Mills gauge theory. The question is what happens to the 
\^ , glueballs in the presence of quark degrees of freedom ? First, a glueball will mix with the ordinary meson with the 

VO ' same quantum numbers so that a pure glueball does not exist in nature. Since the glueball is hidden somewhere in the 

I quark sector, this is one of main reasons why it is so elusive. Second, we shall show that the mass of the pseudoscalar 

ps) ■ glueball could be drastically affected by the dynamic fermion effect. 

P SCALAR GLUEBALL 



It is generally believed that the scalar glueball is hidden somewhere in the isosinglet scalar mesons with masses above 
1 GeV. The argument goes as follows. Many scalar mesons with masses lower than 2 GeV have been observed and 
' they can be classified into two nonets: one nonet with mass below or close to 1 GeV, such as a, K, /o(980) and ao(980) 

that are generally believed to be composed mainly of four quarks and the other nonet with mass above 1 GeV such 
as Kq{143Q), ao(1450) and two isosinglet scalar mesons. This means that not all three isosinglet scalars /o(1710), 
/o(1500), /o(1370) can be accommodated in the qq nonet picture. One of them should be primarily a scalar glueballj^ 
Among the three isoscalar mesons, it has been quite controversial as to which of these is the dominant scalar glueball. 
It was suggested that /o(1500) is primarily a scalar glueball in |[4J, due partly to the fact that /o(1500), discovered in 
pp annihilation at LEAR, has decays to tjtj and rjrj' which are relatively large compared to that of nn and that the 
earlier quenched lattice calculations predict the scalar glueball mass to be around 1550 MeV |5]. Furthermore, because 
of the small production of ;r;r in /o( 17 10) decay compared to that of A'.^, it has been thought that /o( 17 10) is primarily 
ss dominated. In contrast, the smaller production rate of KK relative to ;r;r in /o(1370) decay leads to the conjecture 



' It has been suggested that /o(600) or the a state is a good candidate for the scalar glueball (see e.g. (^). Since a pure glueball state (scalar or 
pseudoscalar) cannot decay into a photon pair or a massless quark pair to the leading order, the identification of /o(600) with a scalar glueball 
seems to be unlikely in view of its broad width, of order 600 ~ 1000 MeV, and its sizable partial width to 77 of order 4keV (^. In the conventional 
tetraquark picture of /b{600), its broad width can be naturally understood as fall-apart decays. 



that /o(1370) is governed by the non-strange light quark content. Based on the above observations, a flavor-mixing 
scheme was proposed [4] to consider the gluebaU and qq mixing in the neutral scalar mesons /o(1710), /o(1500) and 
/o(1370). Best fits to the measured scalar meson masses and their branching fractions of strong decays have been 
performed in several references by Amsler, Close and Kirk 14J, Close and Zhao t^, and He et al. [TJ. A typical mixing 
matrix in this scheme is |@] 

/o(1370)\ /-0.91 -0.07 0.40 \ f\N)\ 

/o(1500) = -0.41 0.35 -0.84 \[\S)\. (1) 
/o(1710)/ \ 0.09 0.93 0.36 / \\G) J 

A common feature of these analyses is that, before mixing, the ss quarkonium mass Ms is larger than the glueball mass 
Mg which, in turn, is larger than the A^(= (mm + dS)l\fl) mass Mat, with Mq close to 1500 MeV and — Mj^ of the 
order of 200 -300 MeV. 

However, the above scenario encounters several difficulties: (i) The isovector scalar meson flo(1450) is now 
confirmed to be the qq meson in the lattice calculation fs*]. As such, the degeneracy of flo(1450) and K^i^WiQ), which 
has a strange quark, cannot be explained if is larger than Mj^ by — 250 MeV. (ii) The most recent quenched lattice 
calculation with improved action and lattice spacings extrapolated to the continuum favors a larger scalar glueball 
mass close to 1700 MeV yj. (iii) If /o( 17 10) is dominated by the content, the decay J ^ (/)/o(1710) is expected 
to have a rate larger than that of J/y/ ^ a)/o(1710). Experimentally, it is other way around: the rate for (b/o(1710) 
production is about 3 times that of J/y/ 0/o(171O). (iv) If /o(1500) has the largest scalar glueball component, 
one expects the T{J/\jf ^ 7/0 (1500)) decay rate to be substantially larger than that of T{J/\jf ^ 7/o(1710)). Again, 
experimentally, the opposite is true. 

In our work [9], we employed two simple and robust results as the input for the mass matrix which is essentially 
the starting point for the mixing model between scalar mesons and the glueball. First of all, we know empirically that 
flavor SU(3) is an approximate symmetry in the scalar meson sector above 1 GeV The near degeneracy of ^^0(1430), 
ao(1470), and /o(1500) has been observed. In the scalar charmed meson sector, D*q(2317) and Dq(2308) have similar 
masses even though the former contains a strange quark. It is most likely that the same phenomenon also holds in the 
scalar bottom meson sector. This unusual behavior is not understood as far as we know and it serves as a challenge to 
the existing hadronic models. Second, an improved quenched lattice calculation of the glueball spectrum at the infinite 
volume and continuum limits based on much larger and finer lattices have been carried out [1]. The mass of the scalar 
gluebaflis calculated tobem(0++) = 1710 ±50 ±80 MeV. This suggests that Mg should be close to 1700 MeV rather 
than 1500 MeV from the earlier lattice calculations f3\. 

We begin by considering exact SU(3) symmetry as a first approximation for the mass matrix. In this case, two of 
the mass eigenstates are identified with ao(1450) and /o(1500) which are degenerate with the mass M. Taking M to 
be the experimental mass of 1474 ± 19 MeV [3], it is a good approximation for the mass of /o(1500) at 1507 ±5 
MeV |3]. Thus, in the limit of exact SU(3) symmetry, /o(1500) is an SU(3) isosinglet octet state and is degenerate 
with ao(1450). In the absence of glueball-quarkonium mixing, /o(1370) becomes a pure SU(3) singlet and /o(1710) 
the pure glueball. When the glueball-quarkonium mixing is turned on, there will be some mixing between the glueball 
and the SU(3)-singlet qq . The mass shift of /o(1370) and /o(1710) due to mixing is only of order 10 MeV. Since the 
SU(3) breaking effect is expected to be weak, it can be treated perturbatively. 



Chiral suppression 

If /o(1710) is primarily a pseudoscalar glueball, it is naively expected that r(G nn)/T[G KK) fa 0.9 
after phase space correction due to the flavor independent coupling of G to PP. However, experimentally there is 
a relatively large suppression of nn production relative to KK in /o(1710) decay, /?(/o(1710)) = r(/o(1710) 
;r;r)/r(/o(1710) KK) = 0.41^Q [y |10] or even smaller To explain the large disparity between Ttn and KK 
production in scalar glueball decays, it was noticed long time ago by Carlson et al. flT'l, by Cornwall and Soni fl2!] and 
revitalized recently by Chanowitz [13] that a pure scalar or pseudoscalar glueball cannot decay into a quark-antiquark 
pair in the chiral limit, i.e., A(G ^ qq) °^ niq. Since the current strange quark mass is an order of magnitude larger than 
niu and m^, decay to KK is largely favored over nn. However, chiral suppression for the ratio Y{G nn) /r(G KK) 
at the hadron level should not be so strong as the current quark mass ratio mu/nis. It has been suggested l\M that niq 
should be interpreted as the scale of chiral symmetry breaking. 

Whether or not G nn is subject to chiral suppression is a controversial issue because of the hadronization 
process from G ^ qq to G ^ nn and the possible competing G qqqq mechanism ifTTl [Til [isi [T6ll . The only 



reliable method for tackling with the nonperturbative effects is lattice QCD. An earlier lattice calculation ITtIi did 
support the chiral-suppression effect with the result r(G ^ rjri) > r(G — > KK) ^ r(G — > TtTt). Although the errors 
are large, the lattice result did show a sizable deviation from the flavor-symmetry limit. If the chiral suppression of 
scalar-glueball decay is confirmed by (quenched or unquenched) lattice calculations, the experimental measurement of 
7?(/o(1500)) = 4.1 ±0.4 alone |3] will be sufficient to rule out the possibility of /o(1500) being a glueball as the ratio 
7?(/o(1500)) should be less than 3/4 due to chiral suppression. Likewise, the identification of /o (600) with a scalar 
glueball is also very unlikely owing to its broad width. 

Guided by the lattice calculations for chiral suppression in G — > PP ijlTll . we have performed a best fit to the 
measured masses and branching fractions. The mixing matrix obtained in our model has the form: 

/o(1370)\ / 0.78 0.51 -0.36\ f\N)\ 

/o(1500) = -0.54 0.84 0.03 |5) . (2) 
/o(1710)/ \ 0.32 0.18 0.93 / \\G) J 

It is evident that /o(1710) is composed primarily of the scalar glueball, /o(1500) is close to an SU(3) octet, and 
/o(1370) consists of an approximate SU(3) singlet with some glueball component 10%). Unlike /o(1370), the 
glueball content of /o(1500) is very tiny because an SU(3) octet does not mix with the scalar glueball. Because 
the nn content is more copious than ss in /o(1710), it is natural that J/\if ^ a)/o(1710) has a rate larger than 
J/xj/ ^/o(1710). Our prediction of T{J/\j/ t(j/o(1710))/r(7/v/ ~^ 0/o(171O)) = 4.1 is consistent with the 
observed value of 3.3 ±1.3 |3]. Moreover, if /o(1710) is composed mainly of the scalar glueball, it will be expected 
that r{J/y/ 7/o(1710)) ^ r(y/v/ 7/o(1500)), a relation borne out by experiment. It is interesting to compare the 
mixing matrices {t^ and In the model of Close and Zhao i^, /o(1710) is dominated by the quarkonium in order 
to explain the suppression of /?(/o(1710)). In our scheme, /o(1710) has the smallest content of ss and the smallness 
of /?(/o(1710)) arises from the chiral suppression of scalar glueball decay. Although /o(1500) in our model has the 
largest content of ss, the KK production is largely suppressed relative to KK, 7? (/o (1500)) = 3[(a/(a + P)]^{Pji/ Pk), 
where /o(1500) — a{\uu) + \dd)) + )3 \ss) and is the cm. momentum of the hadorn h. In SU(3) limit, j3 = —2a and 
this leads to i?(/o(1500)) — 3.9 , in agreement with experiment. 



PSEUDOSCALAR GLUEBALL 

In 1980, Mark II observed a resonance with a mass around 1440 MeV in the radiative J/y/ decay flst ] which was 
subsequently named i(1440) by Mark II and Crystal Ball Collaborations (l^. Shortly after the Mark II experiment, 
i(1440) now known as Tj(1405) was proposed to be a leading candidate for the pseudoscalar glueball. (For an excellent 
review of the £(1420) and i(1440) mesons, see It20il .) Indeed Tj(1405) behaves like a glueball in its productions and 
decays because it has a large production rate in the radiative J/y/ decay and is not seen in 77 reactions. Besides 
77(1405), other states with masses below 2 GeV have also been proposed as the candidates, such as 77(1760) and 
X(1835). 

However, the pseudoscalar glueball interpretation for 77 (1405) is not favored by most of the theoretical calculations. 
For example, quenched lattice gauge calculations predict the mass of the ^ state to be above 2 GeV in |21] and 
around 2.6 GeV in |22, 1]. It is not favored by the sum-rule analysis with predictions higher than 1.8 GeV |23, 24] 
either (for a review of the glueball spectroscopy in various approaches, see Ii25fl ). Hence, we are encountering an 
embarrassing situation that although experimentally 77(1405) is a favored candidate for the pseudoscalar glueball, 
theorists seem to prefer to have a ^ state heavier than the scalar glueball. The motivation of our recent work ll26ll is 
to see if we can learn something about the glueball mass by studying the 77 — 7 7^ — G mixing. 

The 77 — 77' mixing has been well studied by Feldmann, Kroll and Stech 12711 . We extend the FKS formalism to 
include the pseudoscalar glueball G. In the FKS scheme, the conventional singlet-octet basis and the quark-flavor 
basis have been proposed. For the latter, the qq = (uu + dd) / V2 and ss flavor states, labeled by the 77^^ and 77j mesons, 
respectively, are defined. The physical states 77, 77' and G and their decay constants are related to that of the octet, 
singlet, and unmixed glueball states 778, 771 and g, respectively, through the rotation 

1^7) \ / 1^78) \ / f?, ff, \ ( f, f!, \ 

\G) J \\8) ) \ fa fa J fl) 

where = + 54.7° with Q being the 77 — 77' mixing angle in the octet-singlet basis, and is the mixing angle 
between G and 771; that is, we have assumed that 77g does not mix with the glueball. The decay constants for the 



physical and flavor states are defined by 



(o|^7''75?|r?,(^)> = -^fiP""^ my''y5q\ns{p),g{p)) = -^fl,P^ ' 

{Q\sy>'yss\n,{P)) = -if,P^ {0\sY^Y5-'<\ri,{P),8{P)) = -ifl^P" ■ (4) 
Applying the equations of motion for the anomalous Ward identity 

M^frsq) ^2im,qY5q + ^G^yG^\ (5) 

between vacuum and j?]), j?]') and \G), we derive six equations for many unknowns. Hence we have to reply on the 
large Nc counting rules to solve the equations step by step. To the leading order of 1 /Nc expansion, we found that the 
ratio of two of the equations leads to 



ce{s(j) - y^cOAoW^, - sO{c(j) + .JlJlsQAGf-m\ - .JlJlcd^Gm}} y/lf, 
cd{c<j) - yT/3ceAG)fn^,+s9{s<j) - yT/3seAG)^fnl - y^apcml ^ fq ' 



(6) 





where Ag = 1 — cos^g and c(j) (s(j)) is the shorthand notation for cos0 (sin0) and similarly for others. This simple 
equation tells us that the pseudoscalar glueball mass mc can be determined provided that the mixing angle (pa and the 
ratio fs/fq are known. Note that the dependence appears at order of Ac ~ for small ^g. So the solution for mc 
is stable against the most uncertain input ^g. The mixing angles (j) and (pc have been measured by KLOE l,28ilfrom 
the (j) yri,yri' decays. Using the updated results ^ = (40.4 ± 0.6)° and 0g = (20 ±3)° obtained by KLOE1I29I1 in 
conjunction with fs/fq = 1 .352 ± 0.007, we derive the pseudoscalar glueball mass from Eq. ^ to be 

mG = (1.4±0.1)GeV. (7) 

The proximity of the predicted nic to the mass of T](1405) and other properties of Tj(1405) make it a very strong 
candidate for the pseudoscalar glueball. Using the above-mentioned values for (j) and we obtain the rj — rj' — G 
mixing matrix 

0.749 -0.657 0.085 

0.600 0.728 0.331 | | |tji) | . (8) 
-0.279 -0.197 0.940 

Our next task is to check the stability and robustness of our prediction when higher order effects in 1 /Nc are included. 
It turns out that the above simple formula Eq. (|6]l still holds even after keeping the OZI-correcting decay constants, as 
long as they obey the relations fg = V2fg and f^ — f^. Therefore, if excluding the solutions with large and negative 
rri^g, the range ( 1 .4 ± 0. 1 ) GeV of the pseudoscalar glueball mass obtained in Eq. dT) will be more or less respected. 

One may feel very uncomfortable with our solution for niQ as both lattice QCD and QCD sum rules indicate a 
pseudoscalar glueball heavier than the scalar one. The point is that lattice calculations so far were performed under 
the quenched approximation without the fermion determinants. It is believed that dynamical fermions will have a 
significant effect in the pseudoscalar channel, because they raise the singlet would-be-Goldstone boson mass from that 
of the pion to 77 and rj'. Indeed, it has been argued that the pseudoscalar glueball mass in full QCD is substantially 
lower than that in the quenched approximation ll24ll . In view of the fact that the topological susceptibility is large (of 
order 10^^ GeV*) in the quenched approximation, and yet is of order lO^^GeV^ in full QCD and zero in the chiral 
limit, it is conceivable that full QCD has a large effect on the pseudoscalar glueball as it does on rj and rj'. 

Our conclusion of a lighter pseudoscalar glueball is also supported by a recent analysis based on the chiral 
Lagrangian with instanton effects |30]. Two scenarios for the scalar and pseudoscalar glueball mass difference 
AmG = m(0++) -;«(0"+) with the fixed m(0"+) were considered in For 0++ = /o(600) and /o(1710), it is 
found that AmG = -(0.1 --0.3) GeV and (0.2 1.0) GeV, respectively. The first scenario with 0++ = /o(600) cannot 
be realized since there is no any 0^^ glueball candidate at energies between 0.7 1.0 GeV. The second scenario 
indicates the possible candidates of the 0"+ glueball are T7(1295), 77(1405) and Tj(1475). The fact tiiat T7(1405) has 
a large production rate while Tj(1295) has not been seen in the radiative J/y/ decays and that 77(1405) has not been 
observed in 77 reactions while 77 (1475) has supports the proposal that 77(1405) is indeed a good pseudoscalar gluebaU 
candidate. The decay properties of 77(1405) has been recently studied in ll3ll[32ll . 



CONCLUSIONS 



We employed two simple and robust inputs to constrain the mixing matrix of the isoscalar mesons /o(1710), /o(1500), 
/o(1370). In the SU(3) symmetry limit, /o(1500) becomes a pure SU(3) octet and is degenerate with ao(1450), while 
/o(1370) is mainly an SU(3) singlet with a slight mixing with the scalar glueball which is the primary component 
of /o(1710). These features remain essentially unchanged even when SU(3) breaking is taken into account. From the 
analysis of the rj — rj' — G mixing together with the inputs from the recent KLOE experiment, we found a solution 
for the pseudoscalar glueball mass around (1.4±0.1) GeV, suggesting that 77(1405) is indeed a leading candidate for 
the pseudoscalar glueball. Contrary to the mainstream, we thus conjecture that the low-lying pseudoscalar glueball is 
lighter than the scalar one owing to the dynamic fermion or QCD anomaly effect. 

For the lattice community, it will be extremely important to revisit and check the chiral suppression effect in scalar 
glueball decays. The previous lattice calculation in this respect was done almost 15 yeas ago [17.1. If the feature of 
chiral suppression is confirmed by lattice QCD, it will rule out the possibility of /o(1500) and /o(600) as scalar 
glueball candidates. Although it is a difficult and time-assuming task, a full lattice QCD calculation of the ^ 
glueball is desperately needed in order to see if the dynamic fermions will affect the pseudoscalar glueball spectroscopy 
dramatically as they do on tj and Tj'. 
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